Phase properties of operator valued measures in phase space 

T. Subeesh 1,2 and Vivishek Sudhir 3 

1 Dept. of Physics, Indian Institute of Technology Madras, Chennai 600036, India 
2 Dept. of Physics, Amrita School of Engineering, 
Amrita Vishwa Vidyapeetham, Ettimadai 641105, Coimbatore, India^\ 
3 Dept. of Physics, Imperial College London, SW7 2AZ, UK 
(Dated: February 6, 2013) 

Abstract 

The Wigner Phase Operator (WPO) is identified as an operator valued measure (OVM) and 
its eigen states are obtained. An operator satisfying the canonical commutation relation with the 
Wigner phase operator is also constructed and this establishes a Wigner distribution based operator 
formalism for the Wigner Phase Distribution. The operator satisfying the canonical commutation 
relation with the Wigner Phase Operator valued measure (WP-OVM) is found to be not the 
usual number operator. We show a way to overcome the non-positivity problem of the WP-OVM 
by defining a positive OVM by means of a proper filter function, based on the view that phase 
measurements are coarse-grained in phase space, leading to the well known Q-distribution. The 
identification of Q phase operator as a POVM is in good agreement with the earlier observation 
regarding the relation between operational phase measurement schemes and the Q-distribution. 
The Q phase POVM can be dilated in the sense of Gelfand-Naimark, to an operational setting 
of interference at a beam-splitter with another coherent state - this results in a von Neumann 
projector with well-defined phase. 
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I. INTRODUCTION 



The phase of a harmonic oscillator, and hence of the single-mode radiation field, has 
proved to be a challenging question since the inception of quantum theory [IJ |2]. 

The earliest attempts at quantizing the phase, focused on ascribing a consistent meaning 
to the polar decomposition of the field ladder operators [3]. But such a phase operator 
suffers from a pathological non-unitarity which can be traced back to the existence of a 
null-subspace, the one spanned by the ground state, for the number operator jl]. Subse- 
quent attempts relying essentially on the polar decomposition overcomes this difficulty by 
introducing real trigonometric phase operator [SH7] , and phase difference operator [5] . Such 
a phase operator has been found to be useful in an operational setting considered by Noh et 
al. [9]. The hermitian phase operator introduced by Pegg and Barnett [TOj, ITT] , relying on 
well-defined phase states in a finite dimensional Fock space, has since proved to be definitive. 

A complimentary, but intimately related [121 US] approach to the quantum phase is via 
the radial integral of the various quasi-probability distributions in phase space [14, 15J. An 
operator intended to elucidate the phase structure inherent in the Wigner phase distribution 
was constructed and studied in an earlier paper [T6] . 

In [IB] , the operator corresponding to the radially integrated Wigner distribution was 
shown to take the following diagonal form, 



Me) = 




where the state \y/2ue l9 , m) is a particular case of the general state, \z,m) = exp(zat)|m) . 
Then the Wigner phase distribution for the field state p^ is given by, 

Pf{0) = Tr \faM6)) ■ (2) 

Given that the states \z,m) are not orthogonal, but that the Wigner phase operator pw(#) 
obeys the completeness relation, 

f Pw (9)d9 = l, 
Jo 

the proper way to understand this formalism is, as an operator valued measure (OVM) on 
the single-mode field Hilbert space. 

OVM based approach to the quantum phase problem has been considered in the past 
[TTJUH]. The Shapiro-Shepard fTF] approach results in the Susskind-Glogower phase operator 
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as the result of a maximum likelihood estimation on quantum phase measurement, and they 
succeed in showing that the so-constructed OVM has similar measurement statistics as that 
of the Pegg-Barnett hermitian operator. 

In this paper, we find the eigenstates of the Wigner phase OVM (WP-OVM), which hap- 
pen to be rotated quadrature cat states in phase space. Then we construct the operator 
satisfying the boson commutation relation with the WP-OVM, which not surprisingly, is 
seen not to be the usual number operator, in view of the fact that the WP-OVM is not a 
projective operator. The construction of these two operators, satisfying the boson commu- 
tation relation, in the standard oscillator Hilbert space, completes our program of devising 
a Wigner distribution based operator formalism for the Wigner Phase Distribution. 

The latter section is devoted to addressing the problem of non-positivity of the WP-OVM; 
the non-positivity being inherited from the potential negativity of the Wigner function itself. 
We show how a positive OVM can be consistently defined by adopting the view that phase 
measurements are coarse-grained in phase space, i.e., by using a suitable "filter" function to 
mimic the measurement process. For a suitable "filter" function, this leads to a positive OVM 
based on the familiar Q-distribution. This falls in line with earlier observations regarding the 
relation between operational phase measurement schemes and the Q-distribution [T9H22] . 

II. WIGNER PHASE OPERATOR- VALUED MEASURE FORMALISM 

Before embarking on the important task of constructing the eigenspace of the WP-OVM, 
we note a very useful property, viz., 

p w (6)=e^ d p w (0)e-^ e , (3) 

where N = a) a is the usual number operator of the field mode. Hence we see the essential 
connection between rotations on the phase plane and the phase of a given field state; it can 
be shown that this property is also shared by the Pegg-Barnett operator, and can be used to 
establish a weak equivalence between the WP-OVM and the Pegg-Barnett formalism [16J. 
This property also holds for other operators as pointed out by Vaccaro [23] , and so need not 
qualify as a sufficient criterion for a phase operator. Nevertheless, it could be a necessary 
condition. 
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A. Eigenspace of WP-OVM 



Armed with ([3]), we need only consider eigenstates |A) such that, 

p w (0)\X) = X\X), (4) 

so that e |A) are the eignstates of pw{6)- 

Projecting the eigen-equation |4]) onto the position basis (x\, and using the x- 
representation of pwify vi z - ; 

-I /*oo /*00 

p w (0) = — v^ X " dx" / dx' \-x' + x") (x' + x"\, 
2tt J ./_«, 

we obtain the integral equation, 

o 

(x + x) fx(x) dx = An X f\(x') 

,f 

for the ^-representation, f\(x) = (x\X) of the eigenvector |A). Using the technique of taking 
derivatives under an integral, the above integral equation can be cast into the differential 
equation, 

4vrA d 2 J x (x) = f x (-x). 

An elegant way to solve such an equation is to first note that it only connects functions of 
same parity, which suggests expressing the equation in the following form, 

Atx\ (x\p 2 \X) = -(x\V\X), 

where, V — J dx\—x)(x\ is the parity operator, and p = —id x is the usual momentum 
operator. This now implies, 

Since momentum eigenstates do not have definite parity, this equation can only be satisfied 
if | A) are definite-parity superpositions of momentum eigenstates. Careful thought shows 
that the exact form is, 



3/4 



1 \ A 



4tt|A| / . A 



v 



4tt|A| 



(5) 



where the kets are momentum eigenkets, as indicated by the subscript. This actually repre- 
sents two disconnected families of even (for negative A) and odd (for positive A) eigenstates 
of Pvk(O) satisfying, 

/V(0)|±A) = ±A|±A>. 



It is easy to show that these states furnish an orthonormal basis spanning the Hilbert space, 
i.e., 

(A | A') = 8(X - A') 

/oo 
dX |A)(A| = 1. 
-oo 

The Canonical conjugate for the Wigner phase operator, as worked out in Appendix (|A|), 
given by 

Vw(0) = -7r (p 3 x + xp 3 ) V, 

is not compatible with the notion of projective phase measurements and therefore suggests 
that the Wigner phase operator is not a projector. 

III. COARSE-GRAINING AND THE POSITIVE Q PHASE OVM 

Given that the Wigner phase operator is indeed not a true projector, and since it exhibits 
the completeness property of a set of non-orthogonal measurement operators, one is forced 
to view them (taken as a continuum of operators for each value of 9) as an operator valued 
measure. But the serious trouble with such an interpretation is the well known non-positivity 
of the operator [HI [16]. This section is devoted to offering a particular solution to this 
problem, resulting in a positive operator valued measure (POVM). 

The non-positivity of the Wigner phase OVM can be traced back to the non-positivity of 
the postulated definition of the Wigner function for the phase state which has been noted 
previously [HI [TB] . It is well known that the degree of negativity of a Wigner function, the 
so-called non-classicality depth, is operationally equal to the minimum number of thermal 
photons to be added to the state so as to make its Wigner function fully positive [21]. 

Hence, if we take the view that experimentally realized phase measurements are in some 
sense classical, then it could be that one is able to recover a positive operator valued measure 
to describe the situation. So here we suggest that phase measurements are oblivious to the 
connected patches in phase space of area atmost 0(h 2 ) where the Wigner phase distribution 
can be negative; i.e, that classical phase measurements are a coarse-grained sampling [25, 26J. 

That is, we postulate, instead of the phase probability distribution ^ Pp(@) = 
/ W p (£) We(£)d 2 £ with the Wigner function W$ of the phase operator, a coarse-grained 
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phase probability distribution which corresponds to a coarse-grained Wigner function for 
the phase operator, 



where g is the "filter" function in phase space that smears out regions of negativity in Wg. 

Now we need to fix the filter function. To do that, we first note that the phase measure- 
ment scheme, even though postulated to be classical, is assumed to be unbiased, i.e., that 
the measurement process does not introduce any phase-sensitive information in the outcome; 
this implies <?(£') = gd^'D- Next we demand that the phase measurement be optimal, in 
the sense that it does not introduce any additional classical noise. This implies that the 
optimal g is the one which has a support in phase space comparable to the support of a 
minimum uncertainty state, but since it also has to be symmetric, the only possibility is, 
<?(£') = ex P ( — ^2~) ' S enera ^ the fact that a particular form of Gaussian smoothing of the 
Wigner distribution leads to the Q-distribution is very well known [211 123 [28] , and follows 
rather trivially from the Cahill-Glauber continuum of s-ordered distribution functions [29J. 
As a passing remark, we note that this optimal coarse-graining corresponds to the addition 
of a minimum amount of quantum noise as prescribed by the pioneering work of Caves on 
the quantum limits to signal amplification [30] . 

By straightforward computation, we find that such a choice in fact defines a new phase 
probability distribution based on the Husimi Q-distribution, 



so that this can be interpreted as the phase probability distribution arising from the ra- 
dial integration of the Q p . Here, Qg is the angularly peaked Q-distribution satisfying the 
property, 



A. Properties of the Q phase operator 

Following a procedure similar to the one adopted in the construction of the operator 
form of the radially integrated Wigner phase distribution [16], we readily find an operator 






Q p (r cos 9', r sm9')Qg(r cos 6* — 9' , r sin 9 — 9')d9' = Q p (r cos 9, r sin 9). 
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corresponding to the radially integrated Q phase distribution, viz., 

^)^E r h^bo^H (6) 



2vr ^— ' \ 2 y yj n \m\ 

n,m=0 x 

From this expression, we calculate the Q phase probability distribution for Fock states, 

R« .(0) = — , 
l«)H v y 27T 

a uniform distribution, as expected of a meaningful phase operator [TT] . Further, for a 
coherent state we get that, 

where we have defined the rotated quadrature x{6) = 9ft(ae~ ie ); from the properties of the 
error function, it follows that this is a positive phase probability distribution. 
It can also be shown that the Q-phase operator is complete, in the sense that, 

^p Q (9)d9 = l (7) 



o 

In principle, one could evaluate the eigenstates of the Q-phase operator - this leads to a 
Fredholm integral equation with a Gaussian kernel, for the state represented in the coherent 
state basis. The resulting solution appeared to be too tedious and not very instructive, 
therefore we will not present it here. 

As an illustrative application, we consider the phase probability distribution of an even 
cat state, 

M It) + 1-7) 



'2(1 -exp(-2| 7 | 2 )) 

where I7) is the coherent state with complex amplitude 7. It is well known that the Wigner 
phase distribution of this state is in general negative [14J, and thus proves to be difficult to 
interpret classically. In view of the coarse-graining in phase space that was employed herein 
to construct the Q phase operator from its Wigner counterpart, we expect this negativity to 
disappear. Fig. [T] shows the phase space distributions and the associated phase probability 
distributions, clearly showing that the Q-phase probability distribution is positive. 

B. Q phase operator as a POVM 

Very much like the diagonal resolution ([I]), one can arrive at a diagonal form for Pq{0). 
For this case of the Q-phase distribution, this can be most easily deduced by expressing the 




FIG. 1: Figure shows the Q (left panel top) and Wigner (left panel bottom) distributions of the 
even cat state with 7 = 2 (see text for details); regions of negative Wigner function is visible, 
owing to interference in phase space. The right panel shows the Wigner and Q phase probability 
distributions, clearly depicting the lack of negativity in the Q phase probability, in contrast to the 
Wigner phase; inset shows a zoom of the region of interest. 

phase probability distribution as a radial integral over a coherent state expectation of the 
relevant state, i.e., 

1 f 00 

P?(Q) = ~ drr (re i6 \ Pi> \re ie ), 
* Jo 

where \re l6 ) is the phase rotated coherent state exp(iN9)\r) . From this we readily identify, 

1 f°° 

Pq(0) = ~ / drr\re ie }(re w \ } (8) 

which is the required diagonal representation of the operator. 

Together with the completeness result ([7]), and the fact that the Q-phase distribution 
is positive by construction (either as the operator corresponding to the radially integrated 
Q-distribution, or equivalently as a coarse-graining of the Wigner phase operator adopted 
above), this shows that the Q-phase operator (JsJ) is a proper POVM. 

This construction can be directly compared to a previous investigation on the same theme 
[22], where the authors succeed in starting from the radially integrated Q-distribution, and 
then use it to arrive at the Wigner phase distribution and a more general s-parametrized 
phase distribution. In fact, there, and previously [20l [3l]. it is shown that experimentally 
relevant phase measurement schemes with a strong local oscillator are in fact equivalent 
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to the radially integrated Q-distribution, for which we have constructed here, an operator 
representation, as well as a POVM resolution. On the other hand, our construction can 
be compared with a purely operational result obtained in [32], obtained as the result of a 
necessary constraint required of any phase measurement scheme - invariant to phase shifts 
by integer multiples of 2ir. 



C. Dilation of the Q phase POVM 

It is well known, owing to the Gelfand-Naimark theorem, that POVM's on a subsystem 
codify the action of a von-Neumann projector on a bigger system, of which the subsystem 
is a part. This means that a POVM can be understood as the "dilation" of some projector 
into a larger Hilbert space, followed by the action of a unitary operation, succeeded by 
the reduction of the resulting state back into the original Hilbert space. Quite obviously, 
this dilation is ambiguous. But through this ambiguity enters the freedom, and hence the 
generalization, required to specify operational measurement schemes within quantum theory. 

In the case of phase measurements, one plausible operational definition of the measure- 
ment process is the interference of the given state under consideration with a known coherent 
a beam-splitter. In this last part of the paper, we would like to show how the Q 
phase POVM can be understood as the result of a dilation which captures this idea. 

Consider the POVM generated by, 

fl T (/3) = Tr 6 [tj T ® I/?)) ® m Ul] , (9) 

where \(3) is a coherent state in Hilbert space \ip) G T-L a is associated to the projector 
\ip)(ip\, and U T = exp it (atb + dSM is the propagator for a beam-splitter with transmis- 
sivity cos 2 2r. 

(|9p can be explicitly evaluated, assuming a P-representation P(a) for the projector ^(Vl, 
giving, 

r°° r 2n drdcb 

n r (/3) = / / t P{T) 0) |re l< ^ cos t -I- i/3 sin t) (re cos t -\- i/3 sin t| , 

Jo Jo 71 

where the state within the integral is a coherent state of the argument. If this operational 
procedure is to be represented by the Q phase POVM rt8J), then we must have that: (a) 
P(r, (f>) = 5((p — 9), an angularly peaked delta function, and (b) r — > 0, corresponding to a 
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transmissivity of near unity for the state |^>), implying that the ancilla coherent state must 
have p — > oo (for the practical reason of obtaining any signal at all). 

Thus, for an operational setting defined by interference at a beam-splitter with a coherent 
state, the Q phase POVM has the dilation, 

dr d(f> 



PQ{0) 



lim Trb 



TV 



-r5((f)-e)\re i4, ,p)(re i4, ,p\ fj\ 



IV. CONCLUSION 



We have brought to a conclusion, the program of devising an operator formalism for the 
radially integrated phase space distributions. In the process, we discover that the natural 
setting of such an operator formalism is that of an operator valued measure on the single 
mode Hilbert space. 

Although the radially integrated Wigner phase operator is not positive, we show that 
by assuming phase measurements as a phase-unbiased, but quantum-limited, smoothing of 
the Wigner phase operator, one can indeed construct a positive operator valued measure. 
This object is the operator corresponding to the radially integrated Q-distribution, which 
has been previously shown in several contexts to be related to an operational description of 
quantum optical phase. 

Taking the operational approach further, we are able to show that the so-obtained Q 
phase POVM can be dilated into a von Neumann projector in a two-mode Hilbert space. 
By demanding that the operation involved is interference at a beam-splitter with a coherent 
state, we are naturally led to a projector which has a well-defined phase. But quite obviously, 
this is one among an infinite possible dilations of the POVM. 
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Appendix A: Canonical conjugate operator 

Although we have calculated the orthonormal eigenstates of the Wigner phase operator, 
it does not assert that the operator implements projective phase measurements. As has been 
emphasized above, the proper way is to view our operator as an operator valued measure. 
To clearly show the departure of our operator from the phase operators that come under 
the standard von Neumann prescription, we now deduce its hypothetical conjugate operator 
r)w, defined by the canonical commutation, 



If such a conjugate operator existed and were to be compatible with the notion of projective 
phase measurements via pw, then the only possibility, as has been amply pursued in the past 
in the form of the number-phase uncertainty relation 018], is for the conjugate operator to 
be the standard number operator [33] . 

One can explicitly determine f)w; to do so, it is sufficient to find the operator canonically 
conjugate to pw{0)- The latter statement follows by taking into account the manifest am- 
biguity in assigning the orientation of the axes in phase space, and the fact that pw{9) is 
generated from pw(fy by rotations in phase space. 

If |A) are eigenstates of rjw(Q) corresponding to eigenvalue A, then it follows that they 
are given by the delta-normalized wavefunction, 



expression for the eigenstate, one can work out the expression for the conjugate operator, 
which comes out to be, 



[pw,Vw] = i- 




giving the wavefunction in the A-representation, (A|A) 



i 



exp (i A A). By using the above 



Vw(0) 



= —71 



(p 3 x + xp 3 ) V . 
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Since this is not the number operator, it follows that its conjugate pw(Q) does not implement 
projective phase measurements. 
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